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We present results from simulations of the gauge glass model in three dimensions using the parallel 
tempering Monte Carlo technique. Critical fluctuations should not affect the data since we equili- 
brate down to low temperatures, for moderate sizes. Our results are qualitatively consistent with 
earlier work on the three- and four-dimensional Edwards- Anderson Ising spin-glass. We find that 
large scale excitations cost only a finite amount of energy in the thermodynamic limit, and that 
those excitations have a surface whose fractal dimension is less than the space dimension, consistent 
with a scenario proposed by Krzakala and Martin, and Palassini and Young. 
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PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
I. INTRODUCTION 

The ground state structure of spin-glasses is poorly 
understood. While there has, been considerable work on 
Ising-type spin-glass systemsa models with a vector or- 
der parameter symmetry have not yet been analyzed in 
much detail. Here we study the three-dimensional gauge 
glass since it is a simple model with a vector order pararru 
eter in which a finite-temperature spin-glass transitional 
has been well established. 

There are two theories describing the spin-glass phase: 
the "droplet picture" (DP) by Fisher and Husdj and 
replica symmetry breaking picture (RSB) by ParisiE 
According to the droplet picture, a cluster of spins of size 
L costs an energy proportional to L e where 9 is positive. 
This implies that, in the thermodynamic limit, excita- 
tions which flip a finite cluster of spins cost an infinite en- 
ergy. In addition, these excitations have a fractal surface 
with a fractal dimension d s that is smaller than the space 
dimension d. By contrast, RSB follows the exact solution 
of the Sherrington-Kirkpatrick model in predicting that 
there are excitations which turn over a finite fraction of 
the spins and which cost a finite amount of energy in the 
thermodynamic limit. The surface of these excitations 
is space fillingEl, i.e. d s = d. Another difference between 
these models can be quantifi^cLhy looking at the distribu- 
tion of the order parameteiuijiij P(q)- In the droplet pic- 
ture, according to the standard interpretatior£3tL3, P(q) 
is trivial, i.e. there are only two peaks at ±<7EAin the 
thermodynamic limit (^ea the Edwards- Anderson order 
parameter). For finite systems of linear size L, there is 
a tail with weight ~ L~ e down to q = 0. On the con- 
trary, RSB predicts also a tail with a finite weight down 
to q — independent of system size. 

Recently, there have been results biz. Krzakala and 
MartinEl, as well as Palassini and Youngcl (referred to as 
KMPY) for Ising-type systems that find an intermediate 
picture: while large scale excitations cost only a finite 
amount of energy in the thermodynamic limit, their sur- 
face is fractal with d s < d. In this scenario, it is necessary 
to introduce two exponents, 9 and 9' , to describe the sys- 
tem size dependence of the excitation energy, where L e is 
the typical energy of an excitation of size L induced by a 



change in boundary conditions, and 9' describes the size 
dependence of the energy of clusters thermally excited at 
fixed boundary conditions. 

In this paper we test which of the above predictions 
apply to the gauge glass by performing Monte Carlo 
simulations down to low temperatures for a modest 
range of sizes using the parallel tempering Monte Carlo 
methodoiij, as previously done in Ref. |^ for the three- 
and four-dimensional Edwards- Anderson Ising spin-glass. 
We find that P(0) does not decrease with increasing sys- 
tem size and, from data for qi, we deduce that d s < d, 
consistent with the existence of KMPY excitations. 

The layout of the paper is as follows: In Sec. || we 
describe the model as well as the observables measured 
while in Sec. Ill we discuss our equilibration tests for 



the parallel tempering Monte Carlo method. Our results 
are discussed in Sec. [Tv| . In Sec. ^ we summarize our 
conclusions and present some ideas for future work. 



II. MODEL AND OBSERVABLES 



The Hamiltonian of the gauge glass is given by 
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where the sum ranges over nearest neighbors on a square 
lattice in three dimensions of size N — L x L x L and 4>i 
represent the angles of the XY spins. Periodic boundary 
conditions are applied. J is a positive ferromagnetic cou- 
pling between nearest neighbor spins and A^ represents 
the line integral of the vector potential directed from site 
i to site j, 
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and $0 = hc/(2e) is the flux quantum. The Aij are 
quenched random variables uniformly distributed be- 
tween [0,27r]. In this work we set J = 1. Note that 
on average the gauge glass is isotropic even though there 
are local quenched fluxes. 
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This model is often used to describe disordered high-T c 
superconductorstj in a magnetic field since, even though 
it lacks screening, it has the right order parameter sym- 
metry. 

The order parameter of the gauge glass is traditionally 
defined as 

1 N 

i=l 

where we indicate thermal averages by angular brackets, 
(...}, and averages of over the disorder by rectangular 
brackets, [. . .] av . If we average over all possible global 
rotations of the spins then the thermal averages in Eq. (||) 
vanish, so we need to work in an ensemble where global 
rotations of are not permitted (e.g. by applying a small 
symmetry breaking field along the 4> — direction) for 
Eq. (||) to be sensible. This is indicated by the prime on 
the thermal average. In the simulation we evaluate the 
product of the two thermal averages by simulating two 
copies (replicas) of the system with the same quenched 
disorder, and so Eq. (||) becomes 

Q=[(<7>']av, (4) 

where the microscopic spin overlap q is defined by 
1 N 

where a and (3 refer to the two replicas. In the con- 
strained ensemble, if we write q (= q x +iq y ) in terms of its 
real and imaginary parts, then (q y )' = and {q)' = (q x )' . 

In practice it is inconvenient to constrain the ensemble 
by applying a symmetry breaking field, so instead we 
perform a rotation of one replica relative to the other in 
order to maximize q x (and simultaneously set q y to zero). 
It is easy to see that q x in the rotated frame is just |g| 
(which is invariant under rotations). We therefore take 
the spin order parameter to be \q\ and its expectation 
value to be 

Q = [{\q\)\, (6) 

(no prime now indicating an unconstrained thermal av- 
erage) which is to be compared with Eq. (||). In the un- 
constrained ensemble, the probability to get particular 
values for q x , q y only depends on \q\. 

In addition, we will also study the link overlap qi, de- 
fined by 

« = ^ E cos ((# - tf) - ^ - $)) > ( ? ) 

where Nf, = Nd is the number of bonds (d = 3 is 
the space dimension). The sum ranges over all nearest- 
neighbor pairs of spins. Note that while a change in q 
induced by flipping a cluster of spins is proportional to 



the volume of the cluster, qi changes by an amount pro- 
motional to the surface of the cluster. 

The weight in P(\q\) for small |g| varies as L~ s , where 
6' was introduced in section Q. In addition we expect the 
variance of the link overlap to fit to a form Var((ft) ~ 
where, as shown in Ref. || fit = 6' + 2(d — d s ). 

III. EQUILIBRATION 

For the simulatio riS ii w e use the parallel tempering 
Monte Carlo methodtSllj as it allows us to study larger 
systems at lower temperatures. In this technique, one 
simulates several identical replicas of the system at dif- 
ferent temperatures, and, in addition to the usual local 
moves, one performs global moves in which the temper- 
atures of two replicas (with adjacent temperatures) are 
exchanged. This method does not allow us to use the. 
equilibration test first introduced by Bhatt and YoungEZI 
because the system temperature does not stay constant 
throughout the simulation. The equilibration test for 
short range spin-glasses simulated with parallel temper- 
ing Monte Carlo introduced in Ref. ||by Katzgraber et. al 
does not work here either because the disorder is not 
Gaussian. To ensure that the system is equilibrated, we 
therefore require that different moments of q and qi are 
independent of the number of Monte Carlo steps N sweep . 
Figure |l| shows data for several moments of q and qi as a 
function of Monte Carlo Steps. One can clearly see that 
the different moments saturate at the same equilibration 
time. Here we show data for an intermediate size {L = 4) 
since we can better illustrate the procedure by calculat- 
ing longer equilibration times. We also require the ac- 
ceptance ratios of the moves which interchange temper- 
atures to be at least 0.3 or higher and roughly constant 
as a function of temperature. 

As another test for equilibration, we require the dis- 
tribution of q x ,q y , which we call P xy (q x ,q y ), to be sym- 
metric about the origin. Figure || shows a density plot of 
P xy (q x ,q y ) for L — 4 at different temperatures. We see 
clearly that the distributions do not depend on the angle 
from the origin, as required. 

In Table B we show N samp (number of samples), N sweep 
(total number of sweeps performed by each set of spins) , 
and Nt (number of temperature values) , used in the sim- 
ulations. For each size, the largest temperature is 0.947 
and the lowst temperature is 0.050. This is to be com- 
pared withQ'Ea T c « 0.45. The set of temperatures is 
determined by requiring that the acceptance ratios for 
moves which exchange temperatures are satisfactory for 
the largest size, L = 8, and for simplicity the same tem- 
peratures also are used for the smaller sizes. We also 
want the distribution of q at the highest temperature 
to show a Gaussian shape centered at q = to ensure 
that all free-energy barriers have vanished so the sys- 
tem will randomize quickly. The lower-right image for 
T = 0.947 3> T c in Figure || shows this is the case. For 
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all system sizes, the acceptance ratios for global moves 
are always greater than 0.3 for each pair of temperatures. 
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FIG. 1. Moments of the overlap and link overlap, defined 
by Eqs. (^|) and (Q) respectively, as a function of Monte Carlo 
sweeps iVswcop, that each of the replicas perform, averaged 
over the last half of the sweeps. Note that the moments seem 
to equilibrate roughly at the same time and appear to be 
independent of the number of sweeps. The data shown is for 
L = 4 and T = 0.050, the lowest temperature studied. 



L N gamp N gwccp 

3 1.0 x 10 4 6.0 x 10 3 53 

4 1.0 x 10 4 2.0 x 10 4 53 

5 1.0 x 10 4 6.0 x 10 4 53 

6 5.0 x 10 3 2.0 x 10 5 53 
8 2.0 x 10 3 1.2 x 10 6 53 



TABLE I. Parameters of the simulations. AT samp is the 
number of samples, i.e. sets of gauge fields, iV sweep is the total 
number of sweeps simulated for each of the 2Nt replicas for 
a single sample, and Nt is the number of temperatures used 
in the parallel tempering method. 

Since the gauge glass has a vector order parameter 
symmetry, to speed up the simulation we discretize the 
angles of the spins to N$ = 512. This number is large 
enough to avoid any crossover effects to other models as 
discussed by Cieplak et. aO. To ensure a reasonable 
acceptance ratio for single-spin Monte Carlo moves, we 
pick the proposed new angle for a spin within an accep- 
tance window about the current angle, where the size of 
the window is proportional to the temperature T. By 
tuning a numerical prefactor we ensure the acceptance 
ratios for these local moves are not smaller than 0.2 for 
each system size at the lowest temperature simulated. 




FIG. 2. Two-dimensional density plots of Px y (qx,q y ) vs. 
q x , q y for L — 4. The horizontal (vertical) axis represents q x 
(q y ). Dark corresponds to a region of high probability. Note 
that for different temperatures ranging from below to above 
T c the distributions are symmetric around q = 0. (Clockwise 
from top left: T = 0.050, 0.166, 0.520 and 0.947). 
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FIG. 3. Data for the overlap distribution P(|g|) at tem- 
perature T = 0.050 for different system sizes. Note the loga- 
rithmic vertical scale. In this and other similar figures in the 
paper, we only display a subset of all the data points while the 
lines connect all the data points in the set. Thus, the wiggles 
in the lines between neighboring symbols are meaningful. 
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IV. RESULTS 
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FIG. 4. Same as Figure || but at temperature T = 0.199. 

Figures | and | shpw data for P{\q\) at T = 0.050 
and 0.199, well belowO T c « 0.45. There is a clear peak 
for large |g| as well as a tail at small |g|. The weight in 
the tail does not decrease with increasing L, as would be 
expected in the standard interpretations of the droplet 
theory. If anything, the weight actually increases some- 
what for larger sizes. 

Note that P(|<?|) decreases to zero, linearly, at very 
small \q\. This is clearly a phase space factor, since the 
the two-dimensional probability distribution P xy (q x , %), 
plotted in Figure |[ will not diverge for a finite system, 
and P(|<z|) = 2n\q\ P xy (q X7 q y ). In order to have P(|g|) 
tend to a constant for |g| — > 0, a prediction of RSB, the 
region of \q\ over which P(|q|) drops linearly must tend 
to zero for L — * oo, and also P xy (q x , %) must diverge as 
1/ |q in this limit. There is some evidence, particularly 
from Figure||, that P(|g|) stays fiat down to smaller |g| 
for larger L, though the range of sizes is too small to 
make a reliable extrapolation. 

Figures || and p] show data for P(qi). As with the 
distributions of \q\, there is a pronounced peak at large 
^/-values. The width of the distribution decreases with 
increasing system size. 

The variance of P(qi) is shown in Figure ^ for several 
low temperatures. The data is consistent with a power 
law decrease where the (presumably effective) exponent 
varies slightly with T. Extrapolating to T = gives 
(it = 6' + 2{d - d s ) = 0.501 ± 0.04. Assuming 6' « we 
find 



FIG. 5. The distribution of the link overlap at T = 0.050 
for different sizes. Note the logarithmic vertical scale. 




FIG. 6. Same as Figure]^ but at temperature T = 0.199. 



d-d s =0.25 ±0.02, 



(8) 
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implying that system-size excitations have a fractal sur- 
face in the thermodynamic limit as predicted by the 
droplet picture. 
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FIG. 7. Log-log plot of the variance of qi as a function of 
system size L at several temperatures. 



V. CONCLUSIONS 

To conclude, Monte Carlo simulations of the three- 
dimensional gauge glass at low temperatures show that 
the structure of the spin-glass state in this particular 
model agrees qualitatively with previous results for Ising 
spin-glassesu and is in agreement with the KMPY pic- 
ture. From the non-trivial form of P(q) for a mod- 
est range of sizes we infer that system-size excitations 
cost a finite amount of energy in the thermodynamic 
limit. From the variance of the link overlap it ap- 
pears that the surface of these excitations is fractal with 
d — d s = 0.25 ± 0.02. Work is in progress on other vector 
spin-glass models to see if they have the same features 
found here and for the Ising spin-glass. 

These results, however, involve a large extrapolation to 
the thermodynamic limit. There may exist a crossover 
length at larger sizes to a different behavior such as the 
droplet theory or an RSB picture. 
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